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ABSTRACT 

The origin of the magnetic fields in neutron stars, and the physical differences between 
magnetars and strongly magnetised radio pulsars are still under vigorous debate. It 
has been suggested that the properties of the progenitors of neutron stars (the massive 
OB stars), such as rotation, magnetic fields and mass, may play an important role in 
the outcome of core collapse leading to type II SNc. Therefore, knowing the magnetic 
properties of the progenitor OB stars would be an important asset for constraining 
models of stellar evolution leading to the birth of a neutron star. We present here the 
beginning of a broad study with the goal of characterising the magnetic properties of 
main sequence massive OB stars. We report the detection of two new massive magnetic 
stars in the Orion Nebula Cluster: Par 1772 (HD 36982) and NUOri (HD 37061), for 
which the estimated dipole polar strengths, with la error bars, are 1150^200 G an d 



650tf 7 ° Q G respectively. 
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1 INTRODUCTION 

Strong, organised magnetic fields are observed to exist in 
some main sequence stars of spectral type A, B and O. Two 
general models have been proposed to explain the presence 
of these magnetic fields: 

(i) In the dynamo model, the field is generated by a dy- 
namo mechanism, occurring classically in the convective re- 
gions or induced by strong shear during differential rotation. 

(ii) In the fossil model, the field is a remnant from a dy- 
namo active during a previous evolutionary phase, or swept 
up from the interstellar medium (ISM) during star forma- 
tion. This scenario implies that the field must somehow sur- 
vive the various internal structural changes encountered dur- 
ing stellar evolution. The magnetic flux is usually assumed 
to be conserved to some extent. 



* E-mail: Veronique. Petit. l@ulaval.ca 
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Although dynamo models reproduce well the charac- 
teristics of late-type main sequence stars and giants, they 
fail to explain the fields of magnetic early-type stars, as 
their envelopes are primarily radiative. Some models of 
dynamo activity in the small convective cores of those 
stars have been put forward, but they still have fundamen- 
tal difficulties reproducing the observ ed field characteris- 
tics l|Charbonneau fe MacGregorl l2001). Their simple mag- 
netic geometries, lack of significant mass-field strength or 
period-field strength relation, and the fact that the observed 
characteristics of magnetic fields in pre - main sequence Her- 
big Ae/Be stars (Alecian et all 120081; ICatala et all 120071 ; 



iFolsom et aHl2007tlWade et alj|2007tl2005h are qualitatively 
identical to those of their main sequence descendants, point 
toward a fossil origin. Furthermore, the incidence, geome- 
tries and strengths of white dwarf magnetic fields are at 
least qualitatively compatible with evolution from mag- 
netic main sequence A and B stars, suggesting that the 
fields of white dwarfs may also be of fossil origin (e.g. 
IWickra masing he fc Ferrariol 120051 ) . 

In more massive OB stars, magnetic fields have only 
been discovered recently, mostly via clues provided by un- 
usual X-ray properties. Traditionally, the X-ray emission 
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from O and B stars, with a typical level Lx /Lt, i ~ 10 -7 , has 
been explained by r adiative instabilities via a multitude of 
shocks in the wind l|Lucv fc White! Il980l : lOwocki fc Cohen] 
1999). However, the very strong and rotationally mod- 
ulated X-ray emission of the brightest T rapezium star, 
6 1 OriC (Q7, P=15.4d [Cagne et all i|l997l )1 was explained 
bv iBabel fc Montmerld (11997m" in terms of the "magneti- 
cally confined wind shock" model (MWCS). In this model, 
the stellar magnetic field is sufficiently strong, and the ra- 
diative wind sufficiently weak, to allow a dipolar magnetic 
field to confine the outflowing wind in the immediate cir- 
cumstellar environment, resulting in a closed magnetosphere 
with a large-scale equatorial shock which heats the wind 
plasma. In this way, the X-ray emission is enhanced and 
may be modulated by stellar rotation. The MCWS model 
provided a quantitative prediction of a magnetic field in 
8 1 Ori C : such afield (1.1±0 .1 kG) was subsequently discov- 
ered by iDonati et all (l2002h . At the present time, 1 Ori C 
and HD 191612 fl.5kG. IDonati et all [20061 ) are the only 
known O-type stars with directly detected magnetic fields. 
However, it has been speculated that magnetism may be 
widespread among massive stars. Some clues to the presence 
of magnet ic fields comes from X-ray photometry and spec - 
troscopy (|Stelzer et all 120051 ; IWaldron fc Cassinellil 120071) , 
non-themal radio synchrotron emission ( Schnerr et al.l 2007) 
and cyclical variations o f U V wind spectral lines ( Fullertonl 
120031 ; iKaper et ail 1 19961 ). Hence, this lack of magnetic field 
detection may well be due to the fact that direct measure- 
ment of magnetic fields present in the atmosphere of O- 
type and early-B type stars is extremely difficult. These dif- 
ficulties arise from the small number of photospheric lines 
present in the optical spectrum and the large intrinsic width 
of the lines, worsened by the usual fast rotation of these 
stars. 

Neutron stars, evolved from the massive OB-stars, are 
characterised by a wide range of magnetic field strengths. In- 
ferred from spin down rates of radio pulsars, their strengths 
are in the range of 10 11 -10 14 G. Two groups of neutron stars, 
the anomalous X-ray pulsars (AXPs) and the soft gamma 
repeaters (SGRs), host super-strong magnetic fields (10 14 - 
10 15 G), and are referred to as magnetars. It is thought that 
the physical distinction between radio pulsars and magne- 
tars is not simply the dipole field strength, as there is a small 
population of radio pulsar with fields at a magnetar-like 
level, but that does not show the same X-ray characteristics 
(|Kaspi fc McLaughlin! I20051 ). There is some observational 
evidence that neutron stars may evolve from stars as massive 
as 45 Mq , and that many magnetars are linked strongly t o 
these massive stars ( Gacnslcr et al 1 l2005l ; iMuno et al.| [2006) . 

The magnetic flux of 8 1 Ori C (45 M fl ) is (7 ± 3) x 
10 27 G cm 2 (using the stellar radius from ISimon- Diaz et al.l 
(2006)). This magnetic flux is roughly of the same scale as 
the highest field magnetar SGR 18O6-2C0 (~ 3 x 10 28 Gem 2 , 
assuming a 10 km radius) . Therefore, in principle there is 
enough magnetic flux present in a massive magnetic star 
like 8 1 Ori C to explain the super-strong fields seen in some 
neutron stars, under the simple assumption that the mag- 
netic flux is completely conserved during its post-MS stellar 
evolution and transformation into a neutron star. Further- 
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more, provided that OB-star fields are remnants from the 
ISM, the fossil hypothesis could provide a powerful explana- 
tion o f the wide range of magnetic fields present in neutron 
stars |Ferrario fc Wickramas inghe 200iJ). 

On the other hand, it has been suggested that neu- 
tron star magnetic fields could instead be generated during 
the core collapse itself, by a dynamo me chanism induced 
by differential rotation l|Braithwaite! 1200*61 ) . Present studies 
assume that any primordial fields present in the progenitor 
star are weak enough to be expelled by the dynamo process. 
However, if the initial field is strong enough, the evolution 
will be different, as this field is likely to interfere with differ- 
ential rotation and therefore with the dynamo process itself 
jHeger et al.ll2005l ). 

Hence, there seem to be three fundamental parame- 
ters that may play key roles in the origin of neutron star 
magnetic fields, and in the explanation of the differences 
between magnetars and radio pulsars: the primordial mag- 
netic field of the progenitor, the rotation of the star and 
its mass. Therefore, knowing the magnetic properties of the 
progenitor OB stars would be an important asset for con- 
straining models of stellar evolution leading to neutron star 
birth. Many observational efforts are underway to charac- 
terise magnetic fields throughout massive star evolution. We 
present here the beginning of a broad study with the goal 
of characterising the magnetic properties of main sequence 
massive OB stars. 

The Orion Nebula Cluster (ONC) presents a unique op- 
portunity to characterise the magnetic fields of a nearby 
co-evolved and co-environmental population of massive OB 
stars. Furthermore, a Chandra large program, the Chandra 
Orion Ultradeep Pr oject (COUP) was dedicated to observe 
the ONC in X-ray (jStelzer et alj|2005l ). enabling a study of 
the connections between stellar winds, magnetic fields and 
X-rays emission, which will be presented in a subsequent pa- 
per. The ONC contains 9 massive OB stars. They range from 
B3 V (~ 8M ) to 07 V (~ 40 M ), approximately the mass 
range from which neutron stars are thought to be formed. 
In this paper we report the detection of two new massive 
magnetic stars in the Orion Nebula Cluster. 



2 OBSERVATIONS 

We conducted spectropolarimetric observations with the ES- 
PaDOnS spectropolarimeter at CFHT in January 2006 and 
March 2007. We obtained high-resolution (R~65,000) and 
high S/N spectra of 8 of the 9 massive OB stars of the ONC, 
in both epochs. Additional measurements of S 1 OriC and 
Par 1772 were taken with ESPaDOnS in December 2007 and 
with ESPaDOnS's twin Narval, installed at TBL, France, in 
November 2007 respectively. 

A complete circular polarisation observation consists 
of series of 4 sub-exposures between which the polarime- 
ter quarter-wave plate is rotated back and forth between 
position angles, which make it possible to reduce systematic 
errors. For a complete description of observation procedures 
and reduction procedures with the Esprit reduction pack- 
age (which is fundamentally the same as the Libre- Esprit 
package provided by CFHT), see IDonati et al.l (|l997T I. 

In order to increase the magnetic sensitivity of our data, 
we applied the Least Squares Deconvolution (LSD) proce- 
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Figure 1. Least Squares Deconvolved profiles for (^OriC (left), Parl772 (middle) and NU Ori (right). The curves are the mean Stokes I 
profiles (bottom), the mean Stokes V profiles (top) and the N diagnostic null profiles (middle), black for January 2006 and red for March 
2007. 



Table 1. Observation log for the detected stars, along with detection diagnostics and derived longitudinal field components 



Star 


Date (UT) 


HJD 


Total exp. time 


Peak snr a 


LSD snr a 


Detection 


P (%) 


Bj(Q) 


^OriC 


2006-01-09 


53744.792 


4 800s 


1700 


3 000 


Marginal 


99.98 


131 ± 56 


(07V) 


2007-03-10 


54168.835 


3 200s 


1600 


2 600 


Definite 


> 99.99999 


471 ± 53 




2007-12-21 


54456.748 


3 200s 


2 000 


3 700 


None 


66.8 


-53 ± 44 




2007-12-22 


54457.748 


3 200s 


1800 


3 200 


None 


93.2 


122 ± 50 


Par 1772 


2006-01-12 


53747.728 


9 600s 


440 


2 000 


Definite 


99.99997 


-249 ± 77 


(B2V) 


2007-03-07 


54166.699 


9 600s 


760 


3 400 


Definite 


> 99.99999 


84 ±45 




2007-11-11 


54416.550 


6 000s 


370 


1600 


Marginal 


99.93 


-321 ± 95 


NU Ori 


2006-01-12 


53747.852 


9 600s 


1300 


14 000 


None 


50.5 


82 ± 52 


(B0.5V) 


2007-03-08 


54167.703 


9 600s 


1500 


15 000 


Definite 


> 99.99999 


-165 ± 56 



a Per 1.8 km/s pixel for the summed spectra 



dure (|Donati et al.l ll997). which enables the simultaneous 
use of many lines present in a spectrum to detect a mag- 
netic field Stokes V signature. The line masks were carefully 
chosen to exclude Balmer lines and any other lines blended 
with them. Using the LSD technique, we found clear Stokes 
V signatures for 3 stars: the previously-detected 1 Ori C, as 
well as Par 1772 (HD 36982) and NU Ori (HD 37061). The 
LSD profiles for these three stars are shown in Figure [J We 
used the statistical test described bv lDonati et alJ l)l997h to 
diagnose the presence of a signal in either mean Stokes V or 
diagnostic null profiles. A signal is unambiguously detected 
inside the spectral line range whenever the associated detec- 
tion probability is larger than 99.999 per cent (corresponding 
to a false alarm probability smaller than 10 -5 ). The asso- 
ciated longitudinal fi eld can be estima ted with the Stokes 
I and V profiles (e.g. IWade et al.ll2000h . It is important to 
note that the longitudinal field is only given as an indicator 
of the field strength and not the main diagnostic for the pres- 
ence of a magnetic field, because a magnetic configuration 
with a null longitudinal component can still produce a non- 
null Stokes V signature. Table [T] presents an observing log 
for the detected stars, along with the detection probability 
inside the spectral line range and the derived longitudinal 
field component. 

9 1 Ori C is the canonical example of a magnetic O star 
showing rotationally-modulated spectral and X-ray varia- 
tions caused by magnetic confinement of its stellar wind 
l|Babel fc Montmerlel [l997l ; iDonati et aLI l2002t ). This most 
massive star of the Trapezium is a speckle binary with a 
0. 037 separation, composed of a 45 Mq primary (07V) and 
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Figure 2. Longitudinal magnetic field measurements for 6 1 Ori C. 
The M uSiCoS measurements (back circles) are from lWade et all 
(2006) , the UCLE S measurement (red squares) are from 
IDonati et al.l (|2002l ) and the ESPaDOnS measurements (green 
triangles) are from this work. 



a > 6M Q secondary l|Schertl et al.ll2003l ). Optical, UV and 
X-ray features all vary according to a 15.422T0.002 d period 
l|Gagne et al.lll997l ; lstahl et al.lll996l ). We obtained 4 obser- 
vations of this star, and according to this ephemeris, our ob- 
servations correspond to phase 0.49, 0.05, 0.66 and 0.72 for 
January 2006, March 2007 and the 2 observations of Decem- 
ber 2007 respectively. The corresponding derived longitudi- 
nal fields are in good a greement with previous spectropolar i- 
metric measurements ( Donati et al ]|2002l ; IWade et al.ll2006l ) , 
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Figure 3. Marginalized probability densities of the dipole polar 
field strength for Parl772. The magnetic field strength 95% cred- 
ible region is filled in light color. The 68.3% credible region used 
to calculate the lcr error bars is filled in dark color. 



and therefore with their derived dipolar field strength of 
about 1.1 kG (Fig.©. 

Par 1772 is a main sequence (or possibly pre-main se- 
quence) B2 star (~ 6 Mp), with a proje cted rotational ve- 
locity of 80 ± 20km/s (| Wolff et aLll2004h . The March 2007 
Stokes V signature is a good example of a cross-over signa- 
ture, where the longitudinal field component is nearly null 
(here 84 ± 45 G), but the symmetry of the polarized Zee- 
man components is broken by Doppler shifts induced by 
ste llar rotation. The m easurement of 91 ± 193 G obtained 
bv lBagnulo et all (|2006l ) with FORS1 at VLT was likely at 
such a cross-over phase, the magnetic field going undetected 
because of the lower spectral resolving power of that instru- 
ment. 

NU Ori is a triple system, containing a B0.5V primary 
(14 Mo), along with a spectroscopic companion of ~ 3Mq 
(component C, 80 mas separation) and a ~ 1 Mq visual 
companion (componen t B) with a 471 ± 17 mas separation 
jPreibisch et all 1 19991 ). The prima ry is a rapidly r otating 
star, with a wsini = 225±50km/s J Wolff et al.ll2004l ). mak- 
ing it the most rapidly rotating early-B star with a detected 
field. Although such high rotational velocity usually occurs 
only in Be stars, the small and narrow emission in Ha, j3 
and 7 seems more related to nebular emission than to a Be 
behaviour. While there was no formal signal detection for 
the January 2006 observation, the March 2007 observation 
showed a definite signal detection. A close inspection of the 
intensity spectrum revealed the weak, sharp spectral lines 
of the spectroscopic companion. The width of the Stokes V 
signatures compared to the width of the companion spectral 
lines rules out associating the polarization signature with the 
companion - the magnetic field and the profile asymmetries 
are clearly intrinsic to the primary. 



3 SURFACE MAGNETIC FIELDS 

In order to extract the surface field characteristics from the 
observed Stokes V profiles, we compared them with theoret- 
ical profiles for a large grid of dipolar magnetic field configu- 
rations, calculated with the polarised LTE radiative transfer 
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Figure 4. The same as Figure [3] for NU Ori. 

code Zeeman2 (|Landstreetlll988l ; IWade et al.ll200ll ). We sam- 
pled the 4-dimensional parameter space (i, f3, ip, B) which 
describes a centered dipolar magnetic configuration. In such 
a model, i is the projected inclination of the rotation axis to 
the line of sight, /3 is the obliquity of the magnetic axis with 
respect to the rotation axis, tp is the rotational phase and 
B is the polar field strength at the stellar surface. For each 
configuration, we calculated the reduced % 2 of the model fit 
to the observed mean Stokes V profiles. Assuming that only 
the phase may change between two observations of a given 
star, the goodness-of-fit of a given rotation-independent (i, 
P, B)-configuration is expressed in terms of Bayesian prob- 
ability density. This ensures that a good magnetic (i, f3, B)- 
configuration can produce Stokes V profiles that fit all obser- 
vations, as the rotational period is not known with enough 
accuracy to determine a priori the phase difference. Any fea- 
tures that cannot be explained by the rotating dipole model 
are treated formally as additional Gaussian noise, which will 
lead to the most conservative estimates of the parameters, 
according to the maximum entropy principle. 

We can determine the probability density of the field 
strength by marginalizing over inclination and obliquity. We 
then extract a 95% credible region for the surface dipole 
field strength o f each star with the technique described by 
iGregorvl ([2005 ). Figure [3] and Figure U show the resulting 
probability density functions for Par 1772 and NU Ori re- 
spectively. The 95% credible regions are [800, 2450] G for 
Par 1772 and [370, 1220] G for NU Ori. The inferred values 
of the dipole polar strength, with la error bars, are then 
1150+200 G and 650 ±i70 G respectively. 



4 DISCUSSION 

As an illustrative example of the potential of these 
new data to constrain models of neutron star field ori- 
gins, we can compare th e m w ith the predictions of 
iFerrario fc Wickramasinghd (|2006l ) of the magnetic field 
distribution of massive stars (8-45 Mq) on the main se- 
quence. This distribution is based on the observed proper- 
ties of radio pulsars from the 1374-MHz Parkes Multi-Beam 
Survey of isolated radio pulsar, assuming a simple fossil field 
origin with a complete conservation of magnetic flux. They 
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obtained a continuous magnetic field distribution in the pro- 
genitor OB stars, peaking at ~ 46 G with 5 per cenlQ of the 
stars having a field in excess of f kG. 

Of course, our sample contains only 8 stars, but we can 
still make some rough comparisons. Taking the predicted 
field strength distribution, we assume that it is the true 
parent distribution from which we draw a random sample of 
8 stars. We define three possible outcomes: [0-500] G, [500- 
f 000] G and over 1000 G, with respective probabilities de- 
rived from the parent theoretical distribution. According to 
the multinomial distribution, the probability of obtaining 
the distribution of magnetic field strengths observed in the 
ONC is about 1%. 

This result might be interpreted, at first glance, to sug- 
gest that massive OB stars are more magnetic than it would 
be required to explain the magnetic fields of radio pulsars. 
However, some points are important to consider: (i) The 
sample of observed stars may not be representative of a gen- 
eral parent distribution, as the stars all come from the same 
cluster. This region could be unusually magnetic, especially 
if the fields of the OB stars themselves are also of fossil origin 
from the ISM. (ii) The assumed parent distribution is not in 
fact the true parent distribution because some assumptions 
are incorrect, or some elements are missing from the model. 
Examples of such missing physics might be partial flux con- 
servation or the influence of dynamo processes during core 
collapse. 

In order to better explore these possibilities, a larger 
sample of OB stars, from clusters and from the field, must 
be studied in order to increase the population of neutron 
star progenitors with known magnetic properties. Our team 
has undertaken an extensive spectropolarimetric study of 
massive stars in other young star clusters to provide these 
important data. 
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